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THEORY AND OPERATION OF A CARBON DIOXIDE REDUCTION UNIT 


WITH SABATIER REACTION 
By 

Lenwood G. Clark 
ABSTRACT 

A carbon dioxide reduction unit was studied in which carbon dioxide 
and hydrogen are reacted via the Sabatier reaction to produce methane 
and water. The unit is part of the oxygen recovery subsystem of an 
Integrated Life Support System (ILSS) at NASA's Langley Research Center. 
The study included an experimental investigation of the effects of feed- 
gas composition and flow rate, temperature, and pressure on reaction 
efficiency. Consideration is also given to the chemical equilibrium 
characteristics of the Sabatier reaction and an attempt is made to 
correlate experimental reactor performance with theory. 
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V. INTRODUCTION 


Long duration space missions will ultimately require increased 
consideration and use of regenerative life support processes, including 
processes for oxygen recovery. A variety of candidate processes for 
accomplishing oxygen recovery are being developed through the cooperative 
efforts of governmental research laboratories and private industry. Many 
of these processes are surveyed and compared in references 1 to 4. In 
general, oxygen recovery requires two steps: first, removal of the 

metabolic carbon dioxide from the spacecraft's atmosphere and second, 
reduction of this carbon dioxide directly to oxygen or to water which is 
subsequently electrolyzed to produce oxygen. The end waste product of 
most of the physico-chemical processes which accomplish the reduction 
step is elemental carbon, which must be continuously removed from the 
process. Carbon removal problems, however, have been experienced with 
some of these processes, as noted in references 5 and 6, and represent 
a formidable obstacle to the development of reliable life support 
equipment. In addition, potential problems with carbon handling under 
conditions of reduced or zero gravity still remain to be solved . 

One of the most promising carbon dioxide reduction processes is 
based on the "Sabatier” reaction, which does not produce elemental 
carbon. The reaction combines carbon dioxide (COg) and hydrogen (Hg) 
to produce methane (CHj^) and water (HgO) according to the reaction: 

COg + UHg -» 2HgO + CH 4 


1 


2 


Until 10 years ago, this reaction was chiefly considered for use in the 
production of gaseous fuel from coal. Recently, this reaction has been 
considered for use in oxygen recovery processes for space vehicles. As 
envisioned for this application, the product water from the Sabatier 
reaction and some makeup water vould be electrolyzed to produce oxygen 
for crew consumption and hydrogen for the Sabatier reaction. Gas products 
from the Sabatier reaction, vhich would include some unreacted hydrogen 
and carbon dioxide, some uncondensed water, and mostly methane, would be 
vented to space. These losses constitute an important weight trade-off 
consideration and depend primarily on the reaction efficiency. The loss 
of carbon is not an important consideration since it comes from food, 
which would be stored as an expendable for a mission duration suitable 
for using the Sabatier process. 

A variety of catalysts have been investigated, as reported in 
references 7 to 15 , to determine which is most effective in promoting 
the Sabatier reaction to maximum water yield at minimum temperature 
conditions. Both pure metal powder and supported forms of Ruthenium 
and Nickel have been found to be effective catalysts with carbon dioxide 
conversion efficiencies of from 90 to 99 percent reported in the 
literature. Parameters such as temperature, feed-gas composition, 
pressure, reactor length -to-diameter ratio, and feed -gas flow rates 
have been found to affect conversion efficiencies. Several prototype 
oxygen recovery subsystems utilizing this reaction have also been 
developed and investigated to determine operating efficiencies and 
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component interactions . References 14 to 16 report the results of 
these investigations. 

The present investigation was initiated to evaluate the performance 
of the Sabatier carbon dioxide reduction unit which is part of the 
Integrated Life Support System (ILSS) at NASA's Langley Research Center. 
The ILSS is a ground test prototype regenerative life support system 
designed to reclaim and reuse water and oxygen for a four -man crew while 
maintaining a comfortable atmosphere in a test chamber built to simulate 
a spacecraft cabin. A comprehensive summary of the 2-year development 
program is presented in reference 17 with recent test results discussed 
in references 5 and 6 . This study was directed toward determining the 
optimum operating conditions for the ILSS Sabatier reactor, an unusual 
multitubed configuration which uses fluid cooling for temperature 
control . The study included an experimental investigation of the effects 
of feed-gas composition and flow rate, temperature, and pressure on 
reaction efficiency, j i addition, consideration is also given to the 
chemical equilibrium characteristics of the Sabatier reaction and an 
attempt is made to correlate experimental reactor performance with theory. 


VI. LIST OF SYMBOLS 


a,b ,c ,d 


A,B 


C,D 


G 


pO 

G T 


(G$ - Hq) 
T 


% 


K 


M 

P 



R 

R 

S 

T 

T r 

T s 

V 

W 

x 


stoichiometric coefficients 
reactants 

products of reaction 
gas sample 

standard Gibbs' free energy at temperature T, cal/g-mole 
free-energy function, cal/g-mole 

enthalpy of a compound or element at 0° K, cal/g-mole 
chemical equilibrium constant 
molecular weight 

theoretical reaction pressure, psia 
partial pressure of water vapor at the exhaust gas 
saturation temperature, psia 
experimental reaction pressure, psia 
universal gas constant, cal/(g-mole °K) 
volumetric feed-gas ratio, H2 rCOg 
space velocity, hr -1 
theoretical reaction temperature, °K 
experimental reaction temperature, °F 
exhaust gas saturation temperature, °F 
Sabatier reactor core volume (excluding catalyst), ft 5 
weight flow rate of gas or liquid, lb/day 
equilibrium mole fraction 


k 
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x 

X 

AG° 

“?0 

^002 

P 


volume fraction of exhaust gas constituents 
number of moles of C0 2 at equilibrium 
Gibbs' free-energy change in the standard state, 
cal/g-mole 

standard heat of formation of 0° K, cal/g-mole 
carbon dioxide conversion efficiency, percent 

g 

density of gas, lb/ft 


Subscripts 


ch 4 

methane 

O 

o 
r o 

carbon dioxide 

h 2 

hydrogen 

h 2 o 

water 

out 

out of unit 


VII. APPARATUS AMD PROCEDURE 


Test Apparatus 

The experimental investigation was conducted in the NASA-Langley 
Integrated Life Support System (ILSS), which is shown in the photographs 
of figure 1. For this system, the oxygen recovery subsystem was designed 
in accordance with the materials balance shown in figure 2 . Theoretically, 
it is possible to recover 7 .48 lb/day of oxygen (four -man requirement at 
an average of 150-percent basal metabolism rate) from 56 percent of the 
available carbon dioxide and 4.21 lb/day of makeup water. A flow 
schematic of the Sabatier carbon dioxide reduction unit is shown in 
figure 3. For these tests Hg and CO 2 were supplied from high purity gas 
bottles at regulated pressures of about 12 psig. Feed -gas pressures were 
further reduced to about 6 psig by low pressure regulators which 
ordinarily regulate the Hg and COg gas supply pressures from the ILSS 
HgO electrolysis and CO 2 concentration units. Should the delivery 
pressure of either gas fall below 2.0 psig, solenoid valves close 
simultaneously as a safety feature. The desired flow rates of H 2 and 
CO 2 were obtained by adjusting individual flow valves according to the 
pressure drop across calibrated laminar flow elements. The mixed feed 
gas then passes through a desulfurization chamber which removes any 
sulfur or halides present, which could "poison" the Sabatier reaction 
catalyst by forming compounds on the catalyst surface, thereby decreasing 
the number of active reaction sites (refs. 8 , 9 , and l4). Although 
sulfur removal was not necessary for these tests, sulfur contamination 
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in the feed gas is possible from several sources within the ILSS. 

Minute quantities of sulfur generated by the crew could be present in 
the C0 2 delivered to the Sabatier reactor from the molecular sieve unit 
which concentrates metabolic C0 2 from the test chamber atmosphere. In 
addition, one of the 1^0 electrolysis units under investigation uses an 
aqueous sulfuric acid (l^SOi,.) electrolyte. Carryover of this electro- 
lyte through the liquid gas separation membranes into the hydrogen side 
of the electrolysis cells has been a problem as noted in references 5 
and 6. 

Conversion of the Eq and C0 2 to 1^0 and CH^ occurs in the Sabatier 
reactor via the reaction: 

C0 2 + 4Hg -* 21^0 + CH^ + l64o of C0 2 

Some of the exothermic heat of reaction is recovered when the hot exhaust 
products pass back through the desulfurization chamber to warm the feed 
gas . Reaction products as well as unreacted feed gas then pass through 
a condenser where product water is condensed. A porous plate H 2 0 
separator is used to remove water from the unit under conditions of 
reduced or zero gravity. This component was not used, however, and the 
H2O was collected in a trap to permit accurate measurement. Uncondensable 
exhaust products (small amounts of C0 2 , Hg, and K^O, and mostly CH^) are 
vented through a back-pressure regulator to an external exhaust system 
(theoretically to space). 

The desired reaction temperature is maintained by thermostatically 
circulating a silicon fluid (Dow Corning DC-331) at a temperature of 
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375 ° F through the Sabatier reactor. This fluid is supplied at a 
pressure of about 45 psig to the ILSS test chamber by a commercial 
heating and pumping unit and simulates a circuit for delivering waste 
heat to the life support system from a dynamic Brayton cycle radioisotope 
power system. When cooling is required the solenoid valve in the fluid 
circuit is de-energized and closed. The fluid, although at a temperature 
of 375 ° F, is supplied to the reactor as a coolant since reaction 
temperatures are much higher (400° to 700° F). The amount of fluid 
which flows through the reactor can be controlled with a manual valve 
to reduce reactor temperature fluctuations . In addition, an adjustable 
bleed orifice in the solenoid valve bypasses some of the fluid to reduce 
temperature fluctuations of another ILSS subsystem in series on the same 
fluid circuit. When the solenoid valve is closed the available fluid 
flow rate averages 70 lb/hr. When the solenoid valve opens, the fluid 
bypasses the reactor. 

The condenser coolant was an aqueous solution of propylene glycol 
which was held at a nominal temperature of 32° F by a commercial cooling 
and pumping unit. This fluid circuit simulates the low temperature 
coolant fluid circuit which would be used as a heat -transfer loop in 
connection with a space radiator heat sink for temperature control. 

This fluid is supplied at 90 psig with a nominal flow rate of about 
150 lb /hr to the condenser. 

Sabatier reactor and catalyst description .- The Sabatier reactor is 
shown schematically in figure 4 r.nd in the photographs of figure 5 • 

The reactor is basically a cylinder with a length of 12.6 inches and a 
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diameter of 3.72 inches. It is constructed of ASME No. Jl6 stainless 
steel. Smaller in diameter near the bottom, the reactor is sealed with 
a metal conoseal and 3-inch-diameter Marmon flange clamp. Catalyst, in 
tablet form, is contained within the reactor between two perforated 
plates primarily in 12 tubes, which are 0.75 inch in diameter and about 
9 inches long. The perforations in these plates permit gas flow through 
the catalyst bed, but are small enough ( 0.09 inch in diameter) to prevent 
migration of the catalyst. 

The hot DC -331 fluid circulates in the chamber around the catalyst- 
filled tubes. The fluid enters the reactor near the bottom, flows 
upward from side to side around several diverter plates, and exits near 
the top. This fluid is also used for heating purposes while the reaction 
is being initiated. Once initiated, the reaction is exothermic and as 
previously mentioned, cooling is required if the desired reaction 
temperature is less than the passive equilibrium temperature. A sensing 
thermocouple (embedded 3*9 inches into one of the catalyst tubes) and 
optical meter relay control the fluid circuit solenoid valve. Should 
the temperature control mechanism fail, an auxiliary thermoswitch 
(embedded 5*5 inches into one of the catalyst tubes) closes the solenoid 
valve to initiate cooling when the reaction temperature reaches 650° F. 

The 0.0287-cubic-foot volume of the reactor core is filled with 
2.38 pounds (lo 8 l grams) of promoted nickel catalyst. The catalyst is 
Harshaw N-010^ l/ 8 -inch-diameter tablets which come fully activates and 
ready to use as supplied. It is composed of a catalyzing coating 
deposited on a carrier (physical support). The coating is a heterogeneous 
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mixture of about 58-percent nickel (Ni) and 1+2-percent nickel oxide (NiO) 
which represents 80 percent of the entire catalyst mass. The carrier, 
Kieselguhr (diatomaceous earth), comprises the remaining 20 percent of 
the catalyst mass. Although the specific process for preparing N-OIO^ 
is proprietary to The Harshaw Chemical Company, the general method used 
is to deposit a thin coating of nickel oxide on the Kieselguhr and then 
reduce the nickel oxide in a hydrogen atmosphere. The ratio of reduced 
nickel to total nickel for N-OIO^ is 0.60 to O.65. 

Desulfurization chamber and catalyst description .- The desulfuriza- 
tion chamber shown schematically in figure 6 and in the photographs of 
figure 5 is an assembly of concentric 5-0- and 3 . 5 -inch-diameter 
cylindrical tubes constructed of ASME No. }l6 stainless steel. The 
overall unit length is 11.7 inches with sealing at the bottom in the 
same manner as the Sabatier reactor. Desulfurization catalyst, in 
extrusion form, is contained in the inner tube between two perforated 
plates 9*5 inches apart. The hot product gases from the Sabatier 
reactor enter the annular passage between the two cylinders near the 
bottom, circulate around the inner tube, and exit near the top. The 
cool feed-gas mixture of CO2 and H2 flows through the desulfurization 
catalyst bed for sulfur removal and is warmed prior to entering the 
Sabatier reactor. Warming the feed gas acts to reduce the temperature 
gradient along the Sabatier reactor and to reduce the condenser cooling 
requirements . Both the Sabatier reactor and desulfurization chamber 
were insulated with several inches of fiberglass batting. 
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The 0.0371-cubic -foot .volume of the desulfurization chamber core is 
filled with 2.84 pounds (1287 grams) of essentially pure zinc oxide. 

The catalyst is Harshaw Zn-C>4oi 3/l6-inch-long extrusions. This material 
does not actually perform as a true catalyst since in the sulfur removal 
process it gets converted to the corresponding anions of zinc . 

Gas Analysis 

The composition of the inlet and outlet gases to the Sabatier 
reactor was determined by passing samples into a Beckman, Model 320-D, 
process gas chromatograph. Gas analysis is subsequently accomplished by 
comparison of the thermal conductivity of the components of the sample 
gas to the thermal conductivity of a reference carrier gas . This 
chromatograph has four columns for separating the sample into its com- 
ponents . Two identical Carbowax^ 20-M ® columns are to strip out the 
higher molecular weight organic compounds and retard water vapor. Carbon 
dioxide is retained and eluted from a silica gel column which also 
retards water vapor. The last column, a molecular sieve, is used for 
separating hydrogen, oxygen (0 2 ), nitrogen (N 2 ), methane, and carbon 
monoxide (CO). Two samples are taken simultaneously for measurement: 
one 0.3 milliliter and the other 3-0 milliliters in size. The smaller 
sample is used for gases in high concentration (0 to 100 percent by 
volume) and was used for the majority of the measurements during these 
tests. The larger sample is used for the low concentration (0 to 

^Union Carbide trade name. 
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5 percent by volume) gases. A complete analysis of all components In 
both samples required 20 minutes . 

High purity helium (assayed 99*9997 percent by volume) was used for 
the carrier gas at a pressure of 20 pslg at the Instrument ' s metering 
orifice with a nominal flow rate of 36 cubic centimeters per minute. 

The thermal conductivity detector was operated at 100° C with a bridge 
voltage and current of 15 volts and 110 milliamperes, respectively. 
Instrument output was read on a Brown -0.1 to 1.0 millivolt recorder. 

With full-scale deflection equivalent to 100-percent volumetric concen- 
tration of the sample gas constituent, the accuracy of reading the 
recorder chart paper was within ± 0.2 percent. 

The chromatograph was calibrated frequently for all gases being 
analyzed. The primary calibration gas was a mixture: 20.2-percent CO 2 , 

21.1-percent 22.9-percent Ng, 18 . 6 -percent CH 4 > and 17.2-percent CO 
(by volume) . Calibration gases of essentially pure CH^ and H 2 were also 
used. The calibration curves for CO 2 and CHI 4 . (with thermal conductivities 
at 100° C of 5.06 x 10 “5 and 9*99 x 10 cal/(sec-cm° C), respectively) 
are linear and accurate within ±1 percent. Hydrogen, however, with its 
high conductivity of by. 9 k x 10"^ cal/(sec-cm° C) compared to the 
reference gas, helium, with conductivity of 39.85 x 10"^ cal/(sec-cm° C), 
causes poor instrument response. The calibration curve for H 2 is 
nonlinear and measurements vere made with an accuracy of ±k percent. 

Instrumentation 

The desired flow rates of Hg and C0 2 were obtained by adjusting 
manual flow valves according to the pressure drop across calibrated 
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Meriam laminar flow elements. These flow elements were calibrated for 
the particular feed gas at regulated pressure with a Precision Scientific 
wet test meter which is accurate to within ±1 percent. The differential 
pressures for Hq and C0 2 flows were measured with 0 to 5 and 0 to 
10 Inches of water Magnehelic pressure gages with 0.1 and 0.2 inch of 
water divisions, respectively. The reactor pressure was measured with 
a 0 to 5 psig Magnehelic pressure gage with 0.1 psl divisions. 

Iron Constatan thermocouple probes were used for measuring tempera- 
ture. The Sabatier reactor control temperature was read on a 0° to 
750° F Sim-Ply-Trol pyrometer with 10° F divisions. The desulfurization 
chamber temperature (7.5 Inches into the catalyst bed) and the Sabatier 
reactor temperature (5.9 inches into the catalyst bed) were read on a 
75° to 2500° F thermoelectric pyrometer with 50° F divisions. The 
exhaust gas saturation temperature was read on a 0° to 550° F Honeywell 
pyrometer with 5° F divisions . 

Test Procedure 

To verify its integrity prior to testing, the Sabatier C0 2 
reduction unit was pressurized with gaseous N 2 and the leakage rate was 
determined for ambient temperature. This was accomplished by measuring 
the pressure decay rate of a known makeup supply volume while the unit 
pressure waB maintained constant at 10 psig. The total unit leakage 
rate for these conditions was subsequently found to be less than 
5 cc/min. 

Warmup .- Warmup of the Sabatier reactor was accomplished by 
circulating hot DC-351 through the reactor while maintaining a slow 
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C0 2 flow to purge the system of contaminants. When the reaction control 
temperature reached 350° F, the H 2 flow rate was set for about 1.5 lb/day 
and the C0 2 flow rate wa.' adjusted to the corresponding stoichiometric 
value. When the temperature increased to about 350° F, indicating 
reaction initiation, the fluid flow control valve was closed. The 
temperature was then allowed to increase until the desired value was 
reached at which time thermostatic cooling was initiated. Reactor 
warmup to a temperature of 500° F was generally accomplished in 1 hour. 
The selected Hg flow rate was then set by adjusting the appropriate 
valve. An arbitrary CO 2 flow rate was similarly obtained. The desired 
reaction pressure was obtained by adjusting the exhaust gas back-pressure 
regulator. As soon as the reactor control temperature had stabilized, 
test data were obtained. 

Data run .- For a stabilized test condition the inlet C0 2 concentra- 
tion was measured with the gas chromatograph. Assuming the remainder of 
the feed gas to be Hq, the HgiCOjj volumetric feed -gas ratio was then 
calculated to within ±1 percent since the CO? measurement was that 
accurate. The water trap was then drained to a reference mark on the 
sight glass and product water was collected for some time period, usually 
1 hour. During this period the composition of the uncondensable gas 
products being exhausted was measured. Two gas chromatograph readings 
(20 minutes apart) were taken for this purpose and the readings averaged. 
Near the end of the test time period, the Inlet feed-gas ratio was again 
determined to insure that test conditions had not altered. Data runs 
where either the inlet or the outlet readings varied significantly 
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were disregarded. At the conclusion of the data run the product water 
was drained and measured with a graduated cylinder . The C0 2 flow rate 
was then varied and the above procedure repeated until tests for the 
selected condition were completed. 

Data analysis .- The C0 2 conversion efficiency, for a given 

data run can be determined from the expression 


w C o 2 - (w C o 2 ) 


X 100 


where W co and ( W co \ .are the flow rates of C0 2 entering and 
2 V 2 /out 

leaving the system, respectively. The flow rate of the C0 2 feed gas 


entering is 


p C0p W Ho 

* C02 " 7 ~ 


where R is the volumetric feed -gas ratio, H 2 :C0 2 , and Pco 2 > Ph 2 
are the densities of COp and H 2 , respectively, at ambient conditions. 
The flow rate of the C0 2 in the exhaust gas can be calculated from the 
expression 


\ x C0 2 M C0 2 (%g * *C0 2 ~ % 2 o) 

/out x C0 2 M C0 2 + x H 2 M H 2 + x CH 4 M CH 4 f ’fyofatf) 


l6 


where |'» T h - + W C q 2 - is the total flow rate of the exhaust gas 

leaving the system; M^, M^, M CHj+ , M^q are the molecular weights 

of the exhaust gas constituents; and x co 2 ' x H2* x civ x h 2 ° are the 
volume fractions of the exhaust gas constituents. The composition of the 
exhaust gas was determined by the gas chromatograph with the exception of 
the water constituent. The concentration of water can be calculated from 


*H 2 0 

^0 = (P r - 0.3) 


where Pj^q is the partial pressure of the water vapor at the exhaust 
gas saturation temperature, T a , and (P r - 0.3) is the approximate total 
pressure at the system exhaust. 

Another useful design parameter is the space velocity, S, the 
inverse of which is a measure of the dwell time of the reactant gases 
in the reactor core. Space velocity is defined as the hourly volume 
flow of feed gas at standard conditions 1 per unit volume of the reactor 
core (excluding the catalyst) and is given by 


^h 2 


V p h 2 

p co 2i 


Previous investigators (refs. 7 and 8) have found that high space 
velocities, in excess of 500 to 600 per hour, result in reduced C0^ 


Standard conditions will be taken as 70° F, 1^.7 psia which is the 
same as the assumed ambient conditions . 
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conversion efficiency. Space velocities ranging from 174 to 765 per 
hour were investigated during this study. 

Sample data analysis ,- To illustrate the use of the above analysis 
procedure, data recorded during run No. 22 are analyzed. The following 
data were recorded: T r = 500° F, P r = 16 psia, = 1.0 

% 2 0 = * - 5.67, T s = 73° F, Xqq 2 = 0.101, Xjj 2 = 0.046, 

and x CH ^ = 0.812. The densities of Hg and C0 2 at 70° F and 14.7 psia 
are 0.0052 and 0.1144 lb/ft^, respectively. 

The flow rate of C0 2 entering the system for these conditions is 


• _ 0.UIA(1.0) _ lb. 

co 2 0.0052(3.67) " * day 


The partial pressure of the water vapor leaving the system at the 
saturation temperature of 73° F is 0.4019 psia from reference 18. The 
volume fraction of the water constituent is then 


0.4019 „ ^ 

3Ur rt — . - . 0 .026 

^2° (16 - 0.3) 


The flow rate of C0 2 leaving the system is, therefore, 


( W C0 2 ) 


= 0.101(44. 01) (1.0 + 5.99 - 4.24) 

out 0.101(44.01) + 0.046(2. 016) + 0. 8l2(l6. 04) + 0.026(18.02) 


= 0.681 




18 


The C0 2 conversion efficiency for this data run is then calculated 

to be 


1co 2 - * 100 ■ 88 - 6 

which means that of the 5.99 lb/day of C0 2 entering the system, 

88.6 percent was converted to CH^ and . 

The volume of the reactor core is 0.0287 cubic feet and the space 
velocity for this condition is 


s 1_ { 1.0 + _5j99_\ 

24(0.0287) \0. 0052 0.1144/ 


= 555 hr- 1 


and the dwell time is about 10 seconds . 


VIII. ANALYSIS 


Chemical equilibrium exists when the rates of forward and reverse 
reactions are equal (ref. 19). Thus, for a typical reaction 
aA + bB P- cC + dD, there exists an equilibrium constant, K, defined by 


K = 


v c„d 

X C X D 




■jC+d-a-b 


( 1 ) 


where P is the reaction total pressure, x is the equilibrium mole 
fraction, A and B are the reactants, C and D are the products, 
and a, b, c, and d are the stoichiometric coefficients. 

For the Sabatier reaction CO 2 + -» CH^ + 2HqO, the equilibrium 
relation is 


K = 



( 2 ) 


Following the example of reference 8, if the degree of conversion of 
CO 2 is represented by X, then the number of moles of CO 2 at equilibrium 
is 1 - X, of H 2 is 4 - 4x, of CH4 is X, of H 2 0 is 2X, and the total 
present is 5 - 2X. 

For the stoichiometric ratio of 4:1, the equilibrium 

relation becomes 
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X(2X) 2 (5 - 2Xp = X ? (5 - 2X) 2 

(1 - X)(4 - 4x) 4 (5 - 2 X)^P 2 64(1 - X)V 


For the 2.5-percent hydrogen rich condition having a Hg :C0 2 ratio 
of 4.1:1, for example, the equilibrium relation becomes 


X(2X) 2 (5-1 - 2X) 2 
(1 - X) (4 .1 - 4X) 4 ? 2 


From the Gibbs' free-energy relationship enthalpies and standard 
heats of formation, the equilibrium constants may be calculated as a 
function of temperature. This calculation is given in Appendix A for 
reaction temperature from 200° to 800° F and figure 7 shows the result. 

The effects of various parameters on the maximum possible CO 2 
conversion may be determined by relating the mole fraction of COp at 
equilibrium directly to reaction temperature. This can be accomplished 
by calculating the value of the equilibrium constant from the governing 
equation for the selected Hp:C0p ratio for values of equilibrium COp 
mole fraction and reaction pressure. The corresponding reaction 
temperature can then be interpreted from figure 7. 

Figure 8 vas subsequently calculated and shows the maximum possible 
CO 2 conversion at several pressures as a function of reaction temperature 
for a stoichiometric H 2 :CC >2 ratio of 4:1. it can be seen that it is 
desirable to operate the reaction at as low a temperature as possible 
for essentially complete conversion. Low temperature reactions, however, 
are difficult to initiate and sustain. In addition, at low temperatures , 


r 
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condensation of the product water in the reactor could have the 
undesirable effect of decreasing the effectiveness of the catalyst 
(ref. 8). Ideally, it would be desirable to operate the reactor iso- 
themally at 300° F with the inlet gases heated to this temperature. 

In practice, however, there is usually a temperature gradient in the 
reactor bed, with the inlet temperature running higher than the outlet 
temperature because of the exothermic heat release. Conversion would 
subsequently be lower at the reactor inlet and as reactor bed temperature 
decreased at the outlet, conversion would increase. 

Figure 8 also shows that an increase in reaction pressure causes a 
corresponding increase in CO 2 conversion. The effect of pressure becomes 
negligible at temperatures where essentially complete conversion is 
achieved . 

An increase in C0 2 conversion is possible for operation with a 
Hq’.COq ratio greater than 4.0, which would correspond to a hydrogen 
rich condition, as shown in figure 9 • However, this effect becomes 
diminished as the amount of excess hydrogen increases as shown in 


figure 10. 


IX. RESULTS AMD DISCUSSION 


Hie results of experimental studies of the effects of various 
parameters on the COg conversion efficiency of the Sabatier reaction 
are presented in tables 1 and 2. These tables list the various data 
and calculations for 98 of a total 130 data runs, which were made during 
this study. Those data missing were either deleted because of inaccuracy 
or to avoid duplication. Table 1 presents data for 72 data runs at a 
reaction pressure of 16 psia. At this pressure, tests were primarily 
conducted for H2 flow rates of 1.0, 1.5, and 2.0 lb/day with reaction 
temperatures of 4oO°, 500°, and 600° F. The K 2 :C02 feed -gas ratio varied 
from 2.7 to 5.25. No adjustment of gas chromatograph data was made. The 
average total percentage of outlet gas constituents was 100.46 percent 
with the standard deviation being ±1.51 percent, which was well within 
instrument accuracy. Table 2 presents data for 26 data runs conducted 
at a reaction pressure of 20 psia. At this pressure, the H2 flow rate 
was maintained constant at 1.0 lb/day with reaction temperatures of 4oO°, 
500°, and 600° F. The J^:C0 2 feed-gas ratio varied from 2.45 to 4.87. 

The average total percentage of outlet gas constituents was 100 .94 percent 
with the standard deviation being ±1.62 percent, which was again within 
instrument accuracy. 

Inspection of the tabulated data shows that occasionally up to 
15 percent more water is produced than is theoretically possible. It is 
unlikely that this error can be attributed to mistakes in either gas 
analyses , water and time period measurements, or laminar flow element 
calibration. The source of the error was not determined, however, most 
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likely the Hg differential pressure gage was misread and the flow rate 
incorrectly set. If the H 2 flow rate were 15 percent higher than 
anticipated, the calculated CO2 conversion efficiency would be high by 
5 percent. 

The data listed in tables 1 and 2 are presented graphically in 
figures 11 and 12. Some of the data have been replotted in figures 13 
to 15 to better illustrate the effects of selected parameters on C0 2 
conversion efficiency. 


Effect of Excess Hydrogen 

The effect of excess H 2 (above stoichiometric requirements) on C0 2 
conversion efficiency is shown in figures 11 and 12. As predicted by 
theory (fig. 10), conversion efficiency can be seen to vary directly 
with the amount of excess Eq for a constant reaction temperature, with 
the effect becoming diminished as the amount of excess H 2 increases. 

It should be noted that for these figures a small variation in the total 
quantity of gas which passes through the reactor occurs with varying 
excess H2 flow for otherwise constant conditions. Such a small change 
in space velocity, as noted later, does not affect subsequent 
interpretation. 

For a H2 flow rate of 1.0 lb/day, which is the approximate design 
condition, a conversion efficiency approaching 98 percent is possible 
at 600° F for stoichiometric feed conditions at both test pressures as 
shown in figures 11(a) and 12. Conversion efficiency increases slightly 
to about 99 percent for a 3 -percent hydrogen-rich feed condition. 
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At lower reaction temperatures, as noted later, the conversion efficiency 
is lower with an increased amount of excess Hq required to achieve essen- 
tially complete conversion. 

If an amount of excess Hg were available it could be used in one of 
two ways: Either to operate the reactor with hydrogen-rich feed and 

improve conversion efficiency, or to operate the reactor with stoichio- 
metric feed and reduce an additional amount of the available excess COg 
(shown in fig. 2) which would otherwise be vented to space. Obviously, 
the choice lies with the method which produces the largest yield of water. 
The data of figures 11 and 12 tend to indicate that hydrogen-rich feed 
operation of this reactor is not Justified. The data of figure 11(a) 
for a reaction temperature of 500° F can be used to illustrate this 


trade-off as 

follows : 




H 2 :C02 

% 

( lb/ day ) 

*co 2 

(lb/day) 

\o 2 

(percent) 

« 

Wh 2 0 

(lb/day) 

4.0 

1.0 

500 

95 

4.28 

4.4 

1.1 

500 

99 

4.46 

4.0 

1.1 

6.05 

95 

4.70 


If 0.1 pound of excess Hg were available and used to operate the 
reactor with a 10 percent hydrogen-rich feed, 4.46 pounds of water would 
result. However, if this Hp were used to reduce an additional amount of 
CO 2 , 4.70 pounds of water would result, or 0.24 pounds more than pro- 
duced by the other method. However, if all of the H 2 were reacted, 

4.50 and 4.95 pounds of water would result, respectively. It is apparent 
that complete H 2 conversion, and not complete CO 2 conversion, would 
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result in the largest yield of water, since Hg is available in a limited 
q uan tity. The available excess CO2 should therefore be used to operate 
the reactor carbon dioxide-rich to eliminate any unreacted Hg trcm being 
vented to space. 


Effect of Pressure 

The effect of reaction pressure on COg conversion efficiency may be 
determined by comparing figure 11(a) with figure 12. This comparison is 
shown in figure 13 for a Hg flow rate of 1.0 lb/day. Increasing the 
reaction pressure from 16 to 20 psia causes a corresponding increase in 
conversion efficiency with effect becoming diminished as essentially 
complete conversion is achieved. It can be seen that the effect of 
pressure on conversion efficiency agrees with theoretical predictions 
(fig. 8); however, the trend with temperature cannot be correlated, as 
noted later. The effect of pressure is most noticeable in figure 13(a) 
for a reaction temperature of 400° F where for stoichiometric feed 
conditions, conversion efficiency is increased from 91 to 96 percent. 

At 600° F, however, the effect of pressure is negligible with conversion 
efficiency approaching 98 percent for stoichiometric feed conditions. 

At this temperature the reaction pressure could be selected to be con- 
sistent with the pressure requirements of the other units which form the 
IISS oxygen recovery subsystem. 

Effect of Temperature 

The effect of reaction temperature on C0 2 conversion efficiency is 
shown in figure 14. Figures 11(a), 11(c), and 12 also show this effect, 
but less clearly. Comparing figure 14 with figure 9 shows that the 
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effect of temperature with all other factors remaining equal was generally 
as predicted. Chemical equilibrium theory indicates that essenti ally 
conversion is possible up to a given temperature, and that conversion will 
decrease if all of the reactor bed were above this temperature. Figure lU 
shows that increasing the reaction temperature beyond about 600° F causes 
a decrease in conversion efficiency. Operation at a higher temperature 
could produce undesirable product of reaction; namely, carbon monoxide 
and carbon as noted in reference 9* It should be noted that the theory 
is based on an isothermal situation which exists only to a limited degree 
with the present reactor configuration and makes precise temperature 
correlation difficult. A temperature gradient existed along the length 
of the catalyst bed and the reaction temperature , T r , served only as 
a reference measurement. A second thermocouple (embedded 5.9 inches into 
one of the catalyst tubes or 2 inches deeper than the control thermo- 
couple) consistently indicated up to 100° F lower temperature. 

Effect of Flow Rate 

The effect of feed-gas flow rate (expressed as space velocity) on 
COg conversion efficiency is shown in figure 15 for a selected set of 
constant conditions. For space velocities of less than about 600 per 
hour, which would be equivalent to reactor dwell times of greater than 
6 seconds, the effect of flow rate appears negligible. At higher space 
velocities, however, a gradual decrease in conversion efficiency results. 
This is in good agreement with the results of previous investigations 
(refs. 7 and 9 ). 
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Additional Observations 

Several chemical analyses of the product water from the Sabatier 
carbon dioxide reduction unit were made. These analyses showed that 
the water contained up to 45 parts per million (by weight) of nickel 
ions which undoubtedly is the result of the nickel catalyst. Direct 
disassociation of this water in an electrolysis unit to complete the 
oxygen recovery process could present problems over a long period of time. 
If the electrolysis unit vised an acid electrolyte, the nickel ions could 
plate out on the hydrogen electrode. If a basic electrolyte were used, 
the nickel ions could form a flocculant. The effect of contamination 
such as this in the water being processed in an electrolysis unit warrants 
further investigation. The requirement of processing this water in a 
water recovery subsystem for contaminant removal may be necessary prior 
to electrolysis. 

The H£S Sabatier carbon dioxide reduction unit was operated for a 
period of 33^ hours during this investigation. During this period the 
thermostatically controlled fluid circuit which was used for reactor 
temperature control functioned exceedingly well. This temperature con- 
trol feature enabled reaction temperatures to be changed as desired 
during the investigation. However, once optimum operating conditions 
have been determined, a Sabatier reactor could be designed to operate 
at those conditions without the need for such a sophisticated temperature 


control. 
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X. CONCLUDING REMARKS 

An experimental investigation has been made of the effects of various 
parameters on the performance of a Sabatier carbon dioxide reduction unit 
which is part of an Integrated Life Support System. A brief attempt at 
correlating experimental reactor performance with theory was also made. 

The study suggests the following conclusions: 

1. The Sabatier carbon dioxide reduction unit was operated for 
33 6 hours during which time the unit demonstrated the capability of 
operating reliably over a range of feed-gas compositions and flow rates, 
temperatures, and pressures. When the unit was operated with stoichio- 
metric feed, carbon dioxide conversion efficiencies from about 91 to 

98 percent were achieved for "indicated" reaction temperatures from 400° 
to 600° F, respectively. 

2. At all temperatures tested, operation with an excess amount of 
hydrogen feed (above stoichiometric requirements) improved conversion 
efficiency. However, hydrogen-rich operation of a Sabatier reactor as 
part of an oxygen recovery subsystem would not be Justified since com- 
plete conversion of the available H2 would result in the largest yield 
of water. A larger yield of oxygen would subsequently result when the 
water was electrolyzed. 

3. At the optimum operating temperature, 600° F, the effect on 
conversion efficiency of changing reactor pressure from 16 to 20 psia 
was negligible. At lower temperatures where conversion was less com- 
plete, conversion efficiency was higher at the higher pressure. 
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4. The effect of space velocity on conversion efficiency was found 
to be insignificant at space velocities of less than 600 per hour, which 
would be equivalent to reactant gas dwell times in the reactor of greater 
than 6 seconds. 

5. The effects of reaction temperature, pressure, and feed-gas 
composition on carbon dioxide conversion efficiency were found to be 
generally as predicted by chemical equilibrium theory. Since the reactor 
was non-isothermal, however, precise correlation of experimental reactor 
performance with theory was not possible. 
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XIII. APPENDIX A 


DERIVATION OF EQUILIBRIUM CONSTANT, K, 


FROM FREE-ENERGY AND ENTHALPY DATA 


The equilibrium constant, K, is related to the Gibbs' free-energy 
change in the standard state, AG°, according to the formula: 


= R Zn K 


where T is the reaction temperature (°K) and R the universal gas 


constant (1.9870 cal/g-mole °K (ref. 20) ) . The quantity 


may be 


evaluated from 


G§ - Hg 




React 


where 


Grp = standard Gibbs' free energy at temperature T 
= enthalpy of the compound or element at 0° K 


= standard heat of formation at 0° K 


(Grp - Hq) 

Table 5 lists the free-energy function, ^ , and the heat of 

formation, AHS , for the various elements and compounds of interest 
T 0 

for several temperatures. The heat of formation data were obtained 
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from reference 21 and the free-energy function data by linear interpola- 
tion from reference 22. 

Calculation of the equilibrium constant for the Sabatier reaction 
C0 2 + 2H2 -+ 2H2O + CH4. at 477° K (400° F) is as follows: 

= -‘*0.505 - ^222 . 2(40.881 + + 47.251 + ♦ 4(27.596) 

= -40-305 - 33.522 - J21.204 + 47.251 + 196.991 + 110.384 


AG° 

T 


-40.405 


cal 

(g-mole °K) 


, „ Nj° 4q ,4q5 _ 

Zr. R - - jjj 1.^87 20 *355 


-» r 20 • 335 q 
los io K - rw ‘ 8 ' 829 


Therefore, K = 6.745 x 10 8 for 477° K (400° F). 

Table 4- lists the Sabatier reaction equilibrium constants for 
reaction temperatures of 200° to 800° F. These constants were calculated 


in the manner indicated. 


r level 
ripe hatch 
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Figure 1.- The NASA-Langley Integrated Life Suppo: 
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Figure 1.- Concluded- 
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Figure 2.- Oxygen recovery subsystem design mass 
balance with Sabatier mode of CO2 reduction. 




Figure J.- Sabatier COo reduction unit flov schematic. 
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Figure 4.- Sabatier reactor schematic. 
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Figure 8.- Variation of the theoretical CO 2 equilibrium conversion 
with reaction temperature for constant reaction pressures and 
stoichiometric flow. 







TABLE 1.- TABUIATION OF DATA AJD RESULTS FOR 16 PS LA 


48 

















TABLE 1.- Continued 













TABLE 1.- Continued 


50 








TABLE 1.- Continued 


































TABLE 1.- Concluded 


55 










TABLE 2.- TABULATION OF DATA AND RESULTS FOR 20 PSIA 
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constant reaction temperatures . 
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Figure 11.- Concluded, 
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Figure 13.- Effect of excess % on CO2 conversion efficiency for constant 

• lb 

reaction pressures with Wjj 2 = 1.0 


100 



ifo&if&Sl! iil i ii t 4)j ^ »5 ; :i ti**' 










constant H2*.C02 feed ratios with - 1.0 



Effect of flow rate on CO 2 conversion efficiency at 16 psia 
ir stoichiometric at 16 psia flow with T_ = ‘)00° F . 


66 


TABLE 3.- FREE-ENERGY FUNCTION, 
AND HEAT OF FORMATION, AH£ q 


(G? - Hg) 

f » 


, DATA 




Constituents 

T, C 

’K (°f) 

co 2 

% .. 

ch 4 

h 2 o 



(G§ - Hq) 

cal 





T ' 

(g-mole °K) 

366 

(200) 

-45 .071 

-25.757 

-38.061 

-38.730 

394 

(250) 

-45 .694 

-26.305 

-38.719 

-39.371 

422 

(300) 

-46.232 

-26.758 

-39.276 

-39.901 

450 

(350) 

-46.746 

-27.186 

-39.805 

-40 .401 

477 

(400) 

-47.251 

-27.596 

-40.305 

-40.881 

505 

(450) 

-47.742 

-28.013 

-40 .832 

-41.369 

533 

(500) 

-48.183 

-28.364 

-4i .291 

-4l .781 

561 

(550) 

-48.624 

-28.714 

-41.750 

-42.194 

588 

(600) 

-49.061 

-29.050 

-42.178 

-42 .588 

616 

(650) 

-49 .462 

-29.373 

-42.625 

-42 .969 

644 

(700) 

-49.853 

-29.670 

-43.037 

-43.322 

672 

(750) 

-50.243 

-29.968 

-43.448 

-43.674 

700 

(800) 

-50.634 

-30.265 

-43 .860 

-44 .026 




auQ _ 

cal 





(g-mole) 


- 273 

(0) 

-93965 

0 

-15990 

-57107 
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TABLE 4.- SABATIER REACTION EQUILIBRIUM CONSTANTS 


T, °K (°F) 

**10 K 

K 

366 (200) 

14.518 

3.296 X 10 14 

394 (250) 

12.789 

6.151 x 10 12 

422 (300) 

11.296 

1.976 x 10 11 

450 (350) 

9.975 

9.440 x 10 9 

477 (4oo) 

8.829 

6.745 x 10 8 

505 (450) 

7.766 

5.834 x 10 7 

533 (500) 

6.819 

6.591 x 10 6 

561 (550) 

5.956 

9.036 x l(P 

588 (600) 

5.185 

1.531 x 10 5 

616 (650) 

4.466 

2.924 x 10 k 

644 (700) 

3.806 

6.397 x lO^ 

672 (750) 

3.192 

1.556 x 103 

700 (800) 

2.620 

4.093 x 10 2 










